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Abstract. A review of meson emission in heavy ion collisions at incident 
energies around 1 — 2 yl GeV is presented. It is shown how the shape of the 
spectra and the various particle yields vary with system size, with centrality 
and with incident energy. A statistical model assuming thermal and chemical 
equilibrium and exact strangeness conservation (i.e. strangeness conservation per 
collision) explains most of the observed features. 

Emphasis is put onto the study of and K~ emission. In the framework 
of this statistical model it is shown that the experimentally observed equality of 
K'^ and K~ rates at threshold corrected energies \/~s — ^/s^ is due to a crossing 
of two excitation functions. Furthermore, the independence of the to K~ 
ratio on the number of participating nucleons observed between 1 and 10 yl-GeV 
is consistent with this model. The observed flow effects are beyond the scope of 
this model. 



1. Introduction 

Central heavy ion collisions at relativistic incident energies represent an ideal tool to 
study nuclear matter at high densities and at high temperatures. Particle production 
is at all incident energies a key quantity to extract information on the properties of 
nuclear matter under theses extreme conditions. Particles carrying strangeness have 
turned out to be very valuable messengers. At incident energies around 1-2 ^-GeV, 
as available at the SIS accelerator at GSI, Darmstadt, these particles are produced 
below their threshold in NN collisions. Due to strangeness conservation, a meson 
is produced with an associated A, having a threshold of 671 MeV in the centre-of- 
mass system and 1.58 GeV in the laboratory system. K~ is "cheapestly" produced by 
K~K~^ pair creation, yielding thresholds of 987 MeV (cm. system) and 2.5 GeV (lab 
system) . The subthreshold strangeness production is the key parameter governing this 
paper. Furthermore, K~ and behave very differently in nuclear matter. While 
hardly finds a partner to react with, exhibits a strong interaction with nuclear 
matter. 

In this paper, an overview of the shape of the measured meson spectra is given 
and their yields are studied as a function of system size A + A, of the number of 
participating nucleons Apart and as a function of incident energy. These results are 
discussed along with the question whether the observations are in agreement with the 
assumption of a thermal and chemical equilibrium. Of special interest is the ratio 
of / K~ yields as the measured values in heavy ion collisions differ strongly from 
those in pp reactions. These findings have lead to the interpretation that in heavy 
ion collisions the "effective masses" of and K~ might be changed as predicted 
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for dense nuclear matter. It turns out that a statistical approach is very successful to 
explain the measured yields. 

2. Experimental Devices 

At GSI, Darmstadt, three major experimental devices for measuring strange particle 
production are installed: 

• The FOPI setup, a Air detector, designed to study the global features of heavy 
ion collisions at these energies This detector measures all particles and is 
optimized for flow studies and informations related with the reaction plane. It is 
also suited for vertex reconstruction to detect K'^, A, etc. In the present setup 
kaons can only be detected at backward rapidities. A major upgrade is under 
way to improve the detection and trigger capabilities. 

• The magnetic Kaon Spectrometer KaoS specially designed to study rare kaon 
production Q . It is optimized for short flight paths because of the kaon decay. It 
has a very selective trigger to handle high beam intensities. Multiwire chambers 
for tracking are installed to reduce the background in kaon measurements with 

rates 10^ times smaller than for protons. Two detector arrays serve for 
centrality selection and reaction-plane determination. Most measurements are 
performed around midrapidity; in some cases a large fraction of the full phase 
space has been studied by moving the spectrometer. 

• The TAPS detector consists of six arrays of BaF2 modules measuring 7 rays. It 
detects among others the 77 which decays by 7 emission, rj has net strangeness 
zero, but contains ss pairs. 

Since this experimental review is restricted to spectral shapes and yields, mainly 
results from the KaoS collaboration are shown. Very interesting flow studies are given 
in the two following talks by P. Crochet and by Y. Shin. 

3. Spectra and Yields 

Spectra of from mass symmetric systems C+C and Au+Au at incident energies 
from 0.6 to 2.0 A-GeV measured at midrapidity are shown in Fig. ^ exhibiting 
Boltzmann shapes |^ . Their inverse slope parameters increase monotoneously with 
incident energy and the heavier system exhibits harder spectra than the light system 
at the same incident energy. 

Pion spectra show deviations from a Boltzmann shape ^ 0- Their shapes 
are qualitatively understood as composed of two parts, a direct (or thermal) 
pion contribution and pions from decaying baryonic resonances, which produce an 
enhancement at lower kinetic energies |, |, |lO). The slopes of the high -energy 
component agree well with those extracted from the corresponding K'^ spectra. 

Next, we study the inclusive yields. In many cases the phase-space distribution 
has been measured, yielding slight anisotropics. These effects are taken into account 
to obtain the inclusive cross sections and are in the order of at most 30% as compared 
to a mid-rapidity measurement extrapolated to 47r assuming isotropic emission ^. 
These cross sections are divided by the geometrical reaction cross sections (using 
ro=l-2 fm) to obtain multiplicities. For a comparison of different systems these 
multiplicity values are divided by the system mass A as shown in Fig. for 1 and 
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Figure 1. K'^ spec- 
tra measured at midra- 
pidity for Au+Au (up- 
per part) and C+C 
(lower part) at vai|. 
ous incident energies | 
W. The solid lines are 
Boltzmann fits to the 
data. 



1.5 yl-GeV. A striking difference is seen between the trends observed for and 
K~ wliicli exhibit a strong rise of M/A with system size and pions which show a 
shght decrease. The values observed for rj production seem to be rather constant. 
Some ry measurements were not fully inclusive. Those values have been corrected. 
Having in mind that kaons are produced below their NN thresholds, one can think 
of mechanisms causing such a rise: Only multiple collisions can collect the energy 
needed for kaon production and the probability for multistep interactions increases 
with density and/or system size. 

It is of interest to follow how the mass dependence varies with incident energy. 
The upper part of Fig. ^ shows the pion and multiplicities per A for C-l-C and 
Au-|-Au collisions as a function of incident energy. The pion multiplicities {Mt^/A) 
representing the sum of charged and neutral pions Q, are smaller for Au-|-Au than 
for C-l-C whereas the multiplicities per A exhibit the opposite behaviour. These 
trends are already seen in Fig. |2[ With increasing incident energy the differences in 
M/A between the light and the heavy system become smaller as can be seen in Fig. 0. 

The mass dependence of the if + yields is demonstrated in the lower part of Fig.^ 
showing the ratio {M/A)Aa+Au/{M/A)c+G- The enhancement factor for if + reaches 
6 at the lowest measured incident energy. At lower incident energies more and more 
individual NN interactions are needed to accumulate the energy to produce a 
and hence a increasing sensitivity to density occurs. These data represent therefore 
the ideal set of information to extract the nuclear equation of state (EOS). Depending 
on the stiffness of the EOS different densities might be reached. A soft EOS allows 
to compress the nuclear matter easier, higher densities can be achieved and as a 
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Figure 2. Multiplicities per mass number A as a function of A + A for and 
pions {7r+ , tt'', 7r~) (at 1 A GeV) on the r.h.s and (Ml 1 AiGeV) and rj (at 1 
A-GeV) left. Partially preliminary data from [Q, ^ 0, 0- Particles with 

strangeness exhibits a strong rise while those with s=0 do not. 



consequence more can be produced. Yet, qualitative arguments are not sufficient 
to extract the stiffness parameter; transport model calculations are needed and they 
show that these data are in agreement with a soft (fc= 200 MeV) EOS ||rf|. This 
conclusion, however, is put in question by a recent study [f2| . 

Negatively charged kaons are at SIS energies always produced below their 
NN threshold and rather small cross sections are expected in heavy ion collisions. 
Furthermore, K~ can be easily absorbed in nuclear matter. Yet, the measured yields 
turned out to be rather high compared to the yields measured in pp collisions. These 
measurements by the KaoS Collaboration have attracted considerable interest as in 
heavy ion collisions the K~ yield compared to the cross section is much higher 
than expected from NN collisions ||^, |l3|. This is especially evident, if the kaon 
multiplicities are plotted as a function of ^/s — ^/sth where ^/sth is the energy needed 
to produced the respective particle in NN collisions taking into account the mass 
of the associately produced partner. The obvious contrast between NN and AA 
collisions, shown in Fig. ^, has lead to the interpretation of the results by in-medium 
properties which cause e.g. a lower threshold for production when produced in 
dense matter [|7|. Of course, the difference between NN and AA collisions as shown 
in Fig. ^, is not sufficient to conclude on properties of kaons in matter. In heavy ion 
collisions, kaons can be produced by other channels, e.g. ttA K~ N which are not 
available in NN collisions. Only by using detailed transport model calculations one 
can conclude on new properties of kaons in matter [p^ . They demonstrate that by 
including secondary collisions and the newly opened channels the measured K~ yields 
are not reproduced when using bare kaon masses. 
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Figure 3. Upper part: Mul- 
tiplicities of and pions per 
j4 as a function of incident en- 
ergy. The pion data include 
charged and neutral pions. The 
lines represent fits to the data 
(see jij). Lower part: Ra- 
tio of the multiplicities per A 
(Au-|-Au over C-l-C) as a func- 
tion of incident energy. 
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Figure 4. Measured K'^ and 
K~ yields in heavy ion and NN 
collisions as p. function of y/s — 

Vslh ||, |l§ |l4)- < Apart > is 

A/ 2 for heavy ion data and 2 for 
A'^A'' collisions. 



4. Interpretation within a Statistical Model 

The spectral shapes of particles presented so far, are in agreement with the 
interpretation within a statistical concept. The measured multiplicities for pions and 

(Fig. H) evidence a strong contrast; while the pion multiplicity decreases with the 
mass of the colliding system, the K'^ multiplicity rises strongly. The latter observation 
seems to be in conflict with a thermal interpretation, which - in a naive view - should 
give multiplicities per mass number A being constant. 

Usually, the particle number densities or the multiplicities per Apart, here for 
pions, are described in a simplified way by a Boltzmann factor 
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Figure 5. The multiplicity 
of /Apart rises strongly with 
Apart in contrast to the pion mul- 
tiplicity [ pil . This rise can be de- 
scibed by the thermal model in- 
cluding exact strangeness conser- 
vation (see text). 



The production of strange particles has to fulfil strangeness conservation. The attempt 
to describe the measured particle ratios including strange hadrons at AGS and SPS 
using a strangeness chemical potential fj,s is quite successful However, 
this grand-canonical treatement is not sufficient, if the number of produced strange 
particles is small. There, a statistical model has to take care of exact strangeness 
conservation in each reaction as introduced by Hagedorn pO[ |. This is done by taking 
into account that together with each if + another strange particle, e.g. a A is produced: 



with the temperature T, the baryo-chemical potential fj,B, the degeneracy factors 
Qi, the production volume for making the associate pair V (see |2^) and the 
total energies Ei. This formula, simplified for demonstration purpose, neglects other 
combinations leading to the production of K'^ as well as the use of Bose-Fermi 
distributions, which are all included in the computation. The corresponding formula 
for K~ production is similar, but does not depend on /ig. This point will become 
important later on. 



Mk~ ( < Ek- > 
^ exp 



; f Ek+ 

9k+V / -— rrexp 



These formulae lead to a reduction of and yields as compared to the numbers 
calculated without exact strangeness conservation ||2^, Two extreme conditions 
can be seen from Eqs. (2) and (3) ||2^. In the low-number limit, the additional term 
(due to the parameter V) leads to a linear rise of Mx+ /Apart while M^^/ Apart remains 
constant. This is in remarkable agreement with the experimental observations shown 
m Fig. |. For very high temperatures or very large V, the terms in brackets approach 
unity (see Ref. and the formulae coincide with the grand-canonical procedure. 

A further test is the understanding of all particle yields. The measured yields 
(or particle ratios) can be described in this statistical concept by combinations of T 
and hb shown as lines in Fig. ^ for central Ni+Ni collisions at three incident energies. 
Besides the results for ti/ttq all measured particle ratios intersect within a small area 
reflecting T and iib at freeze out. Surprisingly, even the measured K'^ / K~ ratio fits 
into this representation and this ratio does not depend on the choice of the volume 
term V. 
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Figure 6. T vpr siis fig for central Ni+Ni collisions from S yl GeV to 
1.8 j4-GeV. From and new, preliminary data for and K~ |14|. 

As the / ratios are in rather good agreement with the assumption of a 
chemical cquiUbrium, it is of interest to see how the result of the thermal concept, 
using Ref. , appears in the representation as a function of ^/s — ^/sth (as used in 
Fig. which is shown in Fig. 0. Note that the x-axis extends now over a much larger 
range than in Fig. ^. 
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Figure 7. Calcu- 
lated /Apart and 
K~ /Apart ratios in the 
thermal model as a 
function of ^/s — y/sth 
for Ni+Ni collisions. 
The points are results 
for Ni+Ni collisions at 
SIS energies pS, nJ 
and Au+Au at 11X2 
AGeV (AGS) pi]. 
At AGS energies tne 
influence of the system 
mass is negligible. 



At values of — y/sth less than zero, the two excitation function cross. They 
differ at AGS energies by a factor of five. The difference in the rise of the two excitation 
functions can be understood by the formulae given above. The one for production 
contains {E\ — fJ-s) while the other has in the exponent. As these two values are 
different, the excitation functions, i.e. the variation with T, exhibit different rises. 

Furthermore, the two formulae predict that the /K^ ratio for a given collision 
should not vary with centrality as V cancels in the ratio, fndeed, this has been 
observed in Au+Au collisions between 4 and 10.2 A-GeV [|2^, and also at SIS 
energies Two examples are shown in Fig. ^. This independence of centrality 

is rather astonishing as one expects an influence of the different thresholds and the 
density variation with centrality, especially at 1.93 A-GeW where the production 
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Figure 8. Ratio of to K as a function of the number of participants Apa^ 
(left for Au+Au at 10.2 A GeV right for Ni+Ni at 1.5 and 1.93 A GeV Q 
together with the statistical model (dashed line) evidencing the independence of 

*^ Apart 



is above, the K~ production below their respective thresholds. 

Transport-model calculations and other estimates show clearly that strangeness 
equilibration requires a time interval of 40 - 80 fm/c |2^, While the statistical 
model is quite successful describing the particle yields, including strange particles, 
this cannot be taken as proof for chemical equlibration. Especially in case of the K'^ 
production, no strong absorptive channel seems to be available. This point clearly 
deserves further attention. At these low incident energies the strange quarks are 
found only in a few hadrons. The s quark is found only in i^"*", while the s quark 
will be shared between K~ and A (or other hyperons). The latter sharing might be 
in chemical equlibrium as the reactions 

7r° + A p + AT" or 7r~ + A ^ n + 

are strong and have only slightly negative Q-values of -176 MeV. If these reactions 
are the dominating channel, the law of mass action might be applied giving for the 
respective concentrations 

^ - ^. (4) 

As the number of relative to A is small, [A] can be approximated by [K^] and 
rewriting gives 

j|^cxM(^° + vr-)/A (5) 

Figure ^ tests the validity of this relation at SIS and AGS energies showing a 
rather linear rise of K~ /K'^ with M{-k~ + 7r°)/A reflecting a constant value of k. 



5. Beyond the Statistical Model 

While the presented results exhibit an astonishing agreement with the statistical model 
there are clear experimental observations which are beyond the statistical description 
indicating new phenomena. These results are presented by the following talks: 
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Figure 9. The K' /K+ 

ratio as a function of the 
pion multiplicity M{it~ + 
■JT^)/A as a test of the law 
of mass action. Prelimi- 
nary data. 



• The emission of positively charged kaon, measured by the KaoS CoUaboration, 
exhibits an elhptic flow which is hard to understand as the interaction with 
nuclear matter is small Only by including a repulsive K^N potential the 
experimental data can be described. Talk by Y. Shin. 

• The measured pt-dependence of the directed flow of A'+ in Ru+Ru and Ni+Ni 
collisions near threshold, studied by the FOPI collaboration, can only be decribed 
when using a repulsive K^N potential jS^. Talk by P. Crochet. An extremely 
strong antiflow signal of is found in Au+Au collisions at 6 A-GeV by the EOS 
coUaboration Talk by P. Chung. 

Direct experimental evidence for the time evolution of pion emission is presented 
based on the shadowing of spectator matter in certain space-time regions In 
peripheral collisions of Au-fAu at 1.0 A-GeV incident energy the moving spectator 
matter acts like a shutter of a camera shielding the pions, i.e. modifying the pion 
emission pattern according to the spatial distribution of the spectator matter at the 
time of the pion freeze out. The motion of the spectator serves as a calibrated clock. 
This analysis suggests that in An -j- Au collisions at 1 A-GeV most of the high-energy 
pions freeze out within 10 fm/c after the first touch of the nuclei. Low-energy pions 
do freeze out later and over a longer time interval ||33|| . 



6. Summary 

The measured multiplicities of and K~ production at subthreshold energies rises 
more than linearly with system mass A and with the number of participating nucleons 
Apart- The spectra do not show deviations from Boltzmann shapes. 

A global survey of the spectral shapes and yields of the emitted pions, 
and K~ points towards an interpretation within a statistical concept. This model 
takes into account exact strangeness conservation, i.e. the associate production of two 
strange particles (e.g. K'^ and A). It does not only explain the various particle ratios 
(except 77), but also describes the very different Apart-dependences of the K'^ and tt 
multiplicities. 

Despite the apparent success of the statistical model of particle production 
under the assumption of thermal and chemical equilibration and using masses of 
free particles, the present understanding of hadronic interactions contradicts chemical 
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equilibrium for strange particles |Q. While for the channel + N ^ A + n 
is open, for no corresponding channel is available needed to understand how 
equilibrium can be obtained. The kaon yields in the statistical model are based on 
free masses, while transport model largely under-predict the measured K~ yields when 
using bare masses. Only if in-medium masses for K~ are used, which are considerably 
lower, agreement with the measured rates can be achieved. This discrepancy seem 
to put into question our present understanding of interactions at the high densities 
reached in heavy ion collisions. 

Furthermore, the observation of a strong elliptic flow of , and of pronounced pt- 
dependence of a directed flow can only be described using a repulsive K^N potential. 

The discrepancy between the two controversial interpretation needed to describe 
the data might disappear when the collision process is understood as a dynamical 
evolution; in some observables, e.g. particle ratios, the dynamics is lost. However, 
when looking at small areas of the phase space the influence of a dynamical evolution 
could be traced back. 
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